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ABSTRACT 

I argue that, despite repeated claims of Ibata et al., the globular cluster NGC 2419 
does not pose a problem for modified Newtonian dynamics (MOND). I present a new 
polytropic model with a running polytropic index. This model provides an improved 
representation of the radial distribution of surface brightness while maintaining a rea- 
sonable fit to the velocity dispersion profile. Although it may be argued that the 
differences with these observations remain large compared to the reported random 
errors, there are several undetectable systematic effects which render a formal likeli- 
hood analysis irrelevant. I comment generally upon these effects and upon the intrinsic 
limitations of pressure supported objects as tests of gravity. 



Ibata et al. (II) have claimed that observations of the 
distant globular cluster NGC 2419 constitute a severe chal- 
lenge for modfied Newtonian dyanmics (MOND), that this 
object is, in fact, a "crucible" for theories of gravity. This 
claim was based, primarily, upon Newtonian and MONDian 
Michie models in which the phase space distribution is as- 
sumed to be of specific form leading to a density cutoff at a 
finite radius. In my initial response to their paper, I argued 
that these models were too restricted and presented high n 
polytropic MOND models (Sll), in which the radial veloc- 
ity dispersion decreases with radius; to the eye in any case, 
such models appear to be quite consistent with both the 
observed radial distribution of surface brightness and the 
radial profile of velocity dispersion in this cluster. Further, I 
cautioned that it is questionable to rely on formal errors in 
likelihood analyses of observations which may well contain 
observational or intrinsic systematic effects. 

Ibata et al. responded with their own analysis of poly- 
tropes (12) and assert that they had already, in effect, con- 
sidered, and ruled out, a range of non-isothermal MOND 
models which fit the data better than the polytropic model 
that I presented; that, in any case, Newtonian Michie mod- 
els fit the observations significantly better than the MOND 
non-isothermal models. Again this conclusion is based upon 
the use of random errors in a likelihood analysis. It is pri- 
marily this point that I wish to address here because it has 
general relevance to the interpretation of astronomical ob- 
servations. But first, I present a new model in which the 
polytropic index runs with radius. 

12 point to the problem for polytropic MOND models 
as a tension between the predicted surface density distribu- 
tion and the velocity dispersion radial profile. This can be 
summarized, in my own words, as follows: to more precisely 
fit the fall off of surface density within the very small formal 
errors, in particular near the outer cutoff, the polytropic in- 
dex must be rather large (n > 15); that is, the system should 
be closer to isothermal. But to fit the radial decrease in the 
line-of-sight velocity dispersion, the index should be smaller 



(n ~ 10). This result is fairly insensitive to the anisotropy 
radius - the distance beyond which radial orbits begin to 
dominate. 

One obvious solution to this problem is to allow the 
polytropic index to increase with radius (there is nothing 
sacred about the rigid polytropic pressure-density relation) . 
Therefore, I have considered several prescriptions in which 
n increases beyond 100 pc or so, for example, 



n = n exp[r/r p ]. 



(1) 



Below I show the results for one such model. Fig. 1 is the ob- 
served surface density of stars as a function of radius (points) 
and that predicted by the model (solid curve), and Fig. 2 
is the same for the line-of-sight velocity dispersion. For this 
model the central radial velocity dispersion is taken to be 
7.7 km/s, the central density is 38 Mq/pc 2 , the anisotropy 
radius is 11.5 pc, the central polytropic index is 9.2 and the 
scale length for the growth of the polytropic index (r p ) is 850 
pc. Note that with a variable polytropic index, the structure 
equation (Jeans equation) for the run of density contains 
an additional term reflecting the gradient in n. The total 
mass of this model is 6.04 x 10 Mq yielding a mass-to-light 
ratio of 1.4. The value of MOND acceleration parameter is 
ao = 10 -8 cm/s 2 and the MOND interpolation function was 
taken to be standard form as in Sll. 

Now, is this an acceptable match to the observations? It 
is clearly an improvement over my initial polytropic MOND 
model with a fixed n, particularly in matching the surface 
density distribution while maintaining a reasonable repre- 
sentation of the velocity dispersion profile. But Ibata et al. 
will probably tell us, based upon a maximum likelihood anal- 
ysis, that it is completely ruled out by the observations or 
that it is 138 times less probable than the best Newtonian 
Michie model. Such claims, of course, rest on the fact that 
the differences between model and observations, while small 
to the eye, are still larger than the reported measurement 
errors (comparable to the size of the points in the surface 
density distribution). But these statements are misleading 
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Figure 1. Surface brightness distribution in the globular clus- 
ter NGC 2419 (points) compared to the MOND non-isothermal 
(polytropic with running index), non- isotropic model described 
in the text. 
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Figure 2. The observed radial dependence of the line-of-sight 
velocity dispersion in NGC 2419 (points) compared to that of the 
MOND model shown above. 



because they give an impression of precision that is not, 
and cannot be, present in astronomical observations of one 
such object. 

12 point out that in spherical symmetry, given the 
anisotropy factor of a model, /3(r), of a model, there is a 
unique relationship between the run of radial velocity and 
the density distribution. This is certainly true, but the ob- 



served quantities are not the radial velocity or the density 
but the projected versions of these quantities. We know very 
well that small errors in the surface density can translate 
into significant differences in the true density and, hence, in 
the line-of-sight velocity dispersion (12 admit that there are 
problems with de-projection particularly in the outer regions 
where the errors are larger). 

But let us assume that 11,2 perfectly understand the 
systematic and random errors in their observations, and that 
these are vanishingly small. Then there remains a number 
of "known unknowns" - or more precisely "known unknow- 
ables", for example, the symmetry of the cluster. 11,2 argue 
that the cluster appears to be quite round and that this 
supports the assumption of spherical symmetry. But if the 
cluster is oblate or prolate with the symmetry axis lying 
near the line of sight, then the cluster would still appear to 
be round without being spherically symmetric. This would 
clearly alter the relation between the run of line-of-sight ve- 
locity dispersion and the projected density distribution, eas- 
ily by more than the random measurement errors. 

Then there is the unknown of rotation. II consider an 
earlier specific claim of rotation and dismiss this as an effect 
of small number statistics. But they cannot dismiss the pos- 
sibility of rotation in general. Rotation is most likely to be 
present in the outer regions where a systematic velocity of 
about 4 km/s would be sufficient to reduce the model veloc- 
ity dispersion at the outermost measured point in Fig. 2 to 
the observed value. If the rotation axis were within 30 de- 
grees of the line of sight its projection would be less than 2.0 
km/s; given the the lower surface density of stars near the 
cutoff, it is difficult to believe that this would be detectable 
even in the very careful observations of II. 

The point is that conceivable but undetectable intrinsic 
systematic effects vitiate the value of likelihood analyses (I 
haven't even mentioned the "unknown unknowns" - effects 
that we have not thought of ) . If so, it would not the first time 
that sophisticated statistical analyses have been applied to 
astronomical observations dominated by systematic effects. 

There are fewer ambiguities with galaxy rotation curves. 
This is essentially because a two-dimensional object is being 
projected onto a two dimensional sky. There is one compo- 
nent of velocity and the vector most often lies in one plane. 
The optical appearance of the disk combined with a two di- 
mensional radial velocity field can generally unambiguously 
determine the projection, or at least reveal when problems 
- such as non-circular disks, non-planar motion, warping or 
non-circular motion - are present. As is well-known MOND 
works extremely well in predicting the shapes of rotation 
curves using the observed distribution of baryonic matter 
(Sanders & McGaugh 2002) But even so, the rotation curves 
in a few percent of these well-specified systems are not in 
agreement with those predicted by MOND, nor would we 
expect them to be. A perfect theory would not precisely 
predict the rotation curves of all disk galaxies because of 
the intrinsic uncertainties. 

It is worse for spheroidal pressure-supported objects be- 
cause of uncertainties due to projection - shape, rotation - 
as well as the degree and radial dependence of anisotropy 
of the velocity field and, in relaxed systems such as glob- 
ular clusters, the possibility of mass segregation and vary- 
ing mass-to-light ratio . If Ibata et al. had several such ob- 
jects, all showing the same discrepancies, their case would be 
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Comment 3 

stronger, but with only one cluster, the small differences be- 
tween model predictions and observations are hardly defini- 
tive. I point out that opposite conclusion has been reached 
by Scarpa et al. (2011) based upon observations of distant 
globular clusters which appear to be inconsistent with New- 
ton. The point is that strong statements on the nature of 
gravity cannot be made from observations of a single such 
cluster given the uncertainties intrinsic in astronomical ob- 
servations of a three-dimensional object projected onto two 
dimensions. 

I am grateful to Moti Milgrom for helpful comments on 
a draft of this paper, and I thank an anonymous referee of 
my previous paper on this subject for suggesting the idea of 
a running polytropic index. 
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